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Abstract 

In last years a great interest was brought to molecular transport problems 
^ ■ at nanoscales, such as surface diffusion or molecular flows in nano or sub-nano- 

channels. In a series of papers V. D. Borman, S. Y. Krylov, A. V. Prosyanov 
and J. J. M. Beenakker proposed to use kinetic theory in order to analyze 
the mechanisms that determine mobility of molecules in nanoscale channels. 
This approach proved to be remarkably useful to give new insight on these 
issues, such as density dependence of the diffusion coefficient. In this paper we 
revisit these works to derive the kinetic and diffusion models introduced by V. 
D. Borman, S. Y. Krylov, A. V. Prosyanov and J. J. M. Beenakker by using 
classical tools of kinetic theory such as scaling and systematic asymptotic 
analysis. Some results are extended to less restrictive hypothesis. 



to the memory of Carlo Cercignani 
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1 Introduction 



In last years a great interest was brought to micro and nano-flows, in part driven 
by applications such as MEMS, micropumps, lab-on-the-chip systems, carbon nan- 
otubes, molecular sieves, etc. Such flows through micro or nano-geometries (micro- 
or nano-channels or micro or nano-porous materials) exhibit unusual behavior and 
therefore create the need for new or more precise models for numerical simulations. 

The gas flows in micro-geometries are often characterized by rather large Knud- 
sen number but small Reynolds and Mach number. Moreover some phenomena 
which are usually neglected in classical fluid dynamics may take importance, such 
as thermal creep flow, thermal stress slip flow or thermal edge flow. The main tool 
at this scale is the kinetic theory used either for direct numerical simulations or for 
deriving correct fluid limit models able to take into account the main characteristics 
of the flow. A survey of this approach and a complete bibliography can be found 
in ([17]) and ([IS]) and for diffusion models on specific geometries and applications 
the reader can refer to (PQ), ([Sj) arid (0). When considering flows through smaller 
geometries, down to nanoscale, new issues must be addressed (see [3], [H]). The 
first one occurs when the size of the channel is comparable with the range of inter- 
action of the gas molecules with the wall, i.e., for a diameter of a few nanometers. 
Then surface-dominated effects become predominant, and in the vicinity of the sur- 
face, the gas flow loses its three-dimensional character and becomes two- and, in 
the limit, one-dimensional, and the mass flux can dramatically exceeds predictions 
of the Knudsen diffusion model (P3]). Therefore these effects must be included in 
the models. Finally, more complicated effects occur when the size of the pore is 
still smaller (around 0.5 nanometer or smaller), and becomes comparable to the de 
Broglie wavelength and quantum effects must be considered at this scale. 

In this paper we focus on gas flows through geometries with characteristic size 
of some nanometers so that the gas molecule-surface interaction is important but 
the quantum effects of zero-point energy can be neglected. Traditionally, at this 
scale, for the numerical simulation, one uses molecular dynamics (IS]), or for sur- 
face diffusion, lattice gas hopping models. Unfortunately molecular dynamics leads 
to expensive computations and is limited to short space and time simulations. A 
different approach is proposed in a series of papers ( [7], [8], [1], [15], [5], [16], [6]). The 
authors proposed to use kinetic theory in order to analyze the mechanisms that 
determine mobility of molecules in nanoscale channels. This approach proved to be 
remarkably useful to give new insight on these issues and, in particular, it was proved 
that this new kinetic theory of surface diffusion is able to explain in a rather natu- 
ral way the density dependence of the diffusion coefficient which has been observed 
for instance in zeolites. Those articles are not enough known in the community 
of applied mathematics and kinetic theory and it seems interesting to revisit them 
from a mathematical point of view and to derive the kinetic and diffusion models 
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by using classical tools of kinetic theory such as scaling and systematic asymptotic 
analysis. In particular the rigorous approach for modelling molecules trapped by 
surface effects, introduced in (|10j) and ([11]). will be useful in this context. We 
consider here the simple case of molecules moving in a 2D plane. 

The molecule-surface interaction is located on a narrow layer (typically L = 0.3 
nanometers). On microscale geometries, for instance in a pipe with diameter around 
some micrometers the flow of molecules inside this layer need not to be described 
precisely and can simply be modelled by a classical boundary condition such as spec- 
ular or diffuse reflection. For pipes with smaller diameter, about some nanometers, 
but rather long {X = some micrometers or more), the effect of the molecule-surface 
interaction is too much important to be described by a simple boundary condition. 
So the starting point of our study is a kinetic model including the molecule-surface 
interaction through Vlasov terms which takes into account the potential interaction 
of the atoms on the surface and through a relaxation term of molecules by phonons 
which represents the effect of the thermal fluctuations. This is a crude representa- 
tion of the very complex interaction between molecules and the surface but sufficient 
to include in the model much more information than the usual boundary conditions. 
Nevertheless such a model is still very expensive for simulation because of the small- 
ness of the surface layer and the stiffness of the effect of the interaction potential. 
It is therefore useful to look for other models derived from this basic kinetic equa- 
tion by asymptotic analysis. The main influence of the surface on the flow is the 
confinement of some molecules in a narrow zone close to the surface, the "surface 
layer" . On the scale we consider here it is relevant to perform an asymptotic analysis 
when the ratio LjX — > 0. This leads to the second model of the hierarchy, which 
can be seen as a multi-phase model coupling the bulk flow of molecules outside the 
range of surface forces and a two energy group kinetic model describing the surface 
molecules, i.e., the molecules within the range of surface forces (see [7]). 

However the kinetic equations for the surface molecules still contain several time 
scales that can be quite different and make its solution difficult. The first time 
scale is given by tfi the "time of flight" over the potential well, the second one is 
given by r ms the molecule-phonon relaxation time and the last one is given by r 2 , 
the characteristic time for a molecule to cross the surface layer. These times must 
be compared with t max the characteristic time of observation of the overall system. 
When tfi « t max , but r ms « t max , in so far as we are interested in the mass flow 
along the surface we do not need to get information on the variation of the distribu- 
tion function on a short time scale (~ tfi) and it is useful to derive a kinetic model 
at a "mesoscopic" time scale comparable to r ms and t max , by an asymptotic analysis 
when tfi/t max — > 0. If we consider a larger observation time such that r ms « t max , 
the system of surface molecules can be described by a diffusion model obtained by 
an asymptotic limit of the "multi-phase" kinetic model when T ms /t max — > 0. Thus 
various diffusion models are derived according to the respective size of r ms and t z . 
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Kinetic model without molecular collisions 

f(t, x, z, ^ v) 



strong confinment in transverse direction 
(section 2) 



Multi-phase kinetic model 
f(t,x,z,y^v z ) g(t,x,v x , e^ 



large molecule-phonon 
collisions (section 4) 



fast oscillating potential in tangential direction 

(section 3) 



Homogeneized surface kinetic model 
h(t,x,e x ,e z ) 



Surface diffusion models: 
trapped molecules and narrow channel 
N(t,x) 



Figure 1: hierarchy of models. 



These successive rescalings and asymptotic limits lead to a hierarchy of models which 
is sketched on figure EJ 

This paper is organized as follows. The hierarchy of models indicated above is 
derived in section 2 to 4 under a simplifying hypothesis on the surface potential and 
for low density flows for which intermolecular collisions are negligible. In Section 
2 we derive a multiphase kinetic model where the equations for the molecules in 
the range of interaction of the surface reduce from two space-dimension to one- 
space dimension. Section 3 is devoted to the derivation, when tfi/t max — > 0, of a 
mesoscopic kinetic model for the flow of molecules inside the surface layer assuming 
that the surface potential is rapidly oscillating in the x-direction parallel to the 
surface. This homogenized model is derived under the assumption that surface 
molecules can be described by a one energy group of molecules trapped in the surface 
layer. In section 4 we study the diffusion limit,when T ms /t max — > 0, of the kinetic 
model for surface molecules derived in section 2. In a first step we consider that 
surface molecules can be described by a one energy group of molecules trapped in the 
surface layer and we derive the diffusion limit in the isothermal case and afterwards 
in the non-isothermal case. In a second step the diffusion limit is obtained for the 
two energy group model in the configuration of a narrow channel with two surface 
layers and no bulk flow and we consider several regimes according to the ratio of 
r ms and t z . 
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2 Kinetic model for a gas flow with local interac- 
tion with the surface of a solid 



2.1 Introduction 

In this section we address the issue of modeling a flow of molecules at the vicinity 
of a solid wall, for instance in a narrow channel. The influence of the interaction 
of molecules with the surface induces a change in the behavior of the gas flow. It 
has been noticed in previous works ([3], [6]) that the flow loses its three dimensional 
character and becomes two- and, in the limit, one-dimensional and the transport is 
modified by this phenomenon. The same mechanism has been studied in a more 
rigorous way in ([H]). Here, since we assume that the molecules are moving in a 
2D (x, z) plane, the flow near the surface will be a one dimensional flow. Following 
the ideas presented in ([7]) and ([3]), and the mathematical approach proposed in 
(pj]), we set up kinetic models that take into account the effect of the molecule-solid 
interaction and describe the motion of molecules within the range of surface forces 
in lower dimension space. 

For the sake of simplicity we consider that the molecules move in a plane (x, z) 
(see remark 1 for the case of a 3D-flow near a plane wall) and we consider a solid 
occupying the half-space z < 0, and a set of molecules of a gas moving in the 
region z > 0. The state of the gas is described by the distribution function / = 
f(t, x, z, v x , v z ) and its evolution is modeled by the following kinetic equation 

d t f + v x d x f + v z d z f - —d x V d v J - —d z V d v J = I ph + Q m , (1) 

m m 

where V = V(x, z) is the the interaction potential of the molecules with the solid, 
Iph is the molecule-phonon collision integral describing the interaction of the gas 
molecules with the thermal fluctuations of the solid and Q m describes the interaction 
of the gas molecules with each other. 

We assume that the potential V satisfies 

1. < V, 

2. Vx, lim z ^ z) = +oo, 

3. the normal part of the potential is a repulsive-attractive potential, i.e. for 
every fixed x the potential z —> V(x, z) is repulsive (i.e. d z V(x, z) < 0) for 
< z < z m {x) and is attractive (d z V(z) > 0) for z m (x) < z. 

4. The range of the surface forces is finite and thus, the potential satisfies V(x, z) = 
V m for z > L. 
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Figure 2: Repulsive-attractive potential (left) and trajectories of molecules in the 
phase plane (right) 

Following ([15]); we assume a simplified form of the collision term describing 
the interaction between molecules and phonons. This is a crude approximation but 
sufficient to include in the model new physical effects. More precisely we define 



where r ms is the molecule-phonon relaxation time, M is the dimensionless equilib- 
rium distribution function 



m is the mass of a gas molecule, k is the Boltzmann constant and T is the temper- 
ature. Let us remark that the right-hand-side in equation ([!]) satisfies the relation 



which insures the mass conservation. 

Let us notice that the range of the surface forces L is small (typically L pa 0.3nm) 
and in most cases much smaller than the mean free path of the molecules in the gas 
phase. Finally, since the interaction of molecules with phonons is important only 
on a part [0, L*] of the range L of surface forces (0 < z < L*, with L* < L), we can 
assume that I p h = 0, for z > L. 
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Those assumptions imply that in the bulk flow the motion of the gas molecules is 
given by the Boltzmann equation 

d t f + v x d x f + v z d z f = Q rn , (3) 

which is consistent with the classical kinetic theory. 

The main goal in our modeling is to take into account the confinement of some 
gas molecules in the surface layer (0 < z < L) due to the interaction with the surface. 
But since the width of this surface layer is very small, as indicated above, it is useful 
to derive by asymptotic analysis a model by assuming that the ratio L/X tends 
to zero (where X is the characteristic length of observation). The phenomenon of 
confinement is mainly transverse to the boundary, but since the interaction potential 
works in both x and ^-directions the analysis is rather technically complicated. We 
suggest to simplify once more the configuration to make easier the analysis while 
keeping in the model enough important features to put in evidence new phenomena 
near the surface at the nanoscale. So we assume that the motion of the gas molecules 
inside the bulk flow is given by (E]) and that the motion of the gas molecule inside 
the range of the surface forces can be separated into two independent components, 
the tangential and the transverse ones. This can be achieved by assuming that 

d x V(x, z) = f(x) and d z V(x, z) = i/j(z), 

which can be obtained by assuming that the interaction potential writes, for z < L 

V(x,z) =U(x) + W(z), (4) 

where W is a repulsive attractive potential with W(z) = W m for z > L and where 
U, the tangential component of the potential, satisfies < U{x) < U m . Then if 
W m » U m , we can assume that 

for z>L, V{x,z) = U{x) + W m nW m , (5) 

which is consistent with the assumption that the range of the surface forces is smaller 
than L (see figure fl2])). 

Another simplification is related to the intermolecular collisions. The study of 
the influence of the surface layer on the flow is relevant only in geometries with one 
characteristic dimension (for instance the size of the pore in a porous medium or 
the diameter of a pipe) about some nanometers. Thus in a first step it is reasonable 
to assume that we are working at a scale where the intermolecular collisions can be 
neglected so that 

Qm = 0. (6) 

These assumptions (HUME]) are used to derive a hierarchy of models describing a 
gas flow in the vicinity of the surface. 
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Remark 1. It should be noted that, under the assumptions (0|) and (0|) ; we can 
assume that the molecules are moving in the original 3D space (x, y, z)(z > 0). More 
specifically, if we assume the 3D version of |IP originally for the molecule-phonon 
interaction, then f corresponds to the marginal velocity distribution function, i.e., 
the integral of the original velocity distribution function with respect to v y from — oo 
to oo. 



2.2 Motion of molecules in a repulsive-attractive potential 

Before all we need to introduce some notations associated with the motion of a 
molecule in a one- dimensional repulsive-attractive potential. Let us denote by e z 
the total energy of the normal motion 

e z = -mv 2 z + W(z), 
and let us introduce the "equivalent velocity" 



for v z ^ 0, e z = e z (z,v z ) = sgn(v z )\jv 2 z + —W(z), 

m 



so that e z is an odd function of v z and 



lim Vz _>Q + e z {z,v z ) = \I—W(z) ^ lim Vz ^ Q - e z (z,v z ) = -J—W(z). 

m \ m 

We denote z m , < z m < L, the unique point where W(z m ) = (see figure ([2])), and 
we define Z-(e z ) and z + (e z ) in the following way: 

• for e z 7^ 



< z-(e z ) < z m < z + (e z ), 
for < y/2W m /m < \e z \, z + (e z ) = L, 

W{z_{e z )) = -el 

for < \e z \ < y/2W m /m, W(z + (e z )) = ^e 2 z , 

W(z.(e z )) = ^el 

• for e z = 0, 

z_(0) = z m = z + (0) 

We notice that the particles with equivalent velocity \e z \ < y^2W m /m are trapped 
(i.e. cannot leave the "surface layer" < z < L) and move between z_(e z ) and 
z + (e z ), but the particles with \e z \ > ^2W m /m are free and move between z_(e z ) 
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and z + (e z ) = L and can go out the surface layer and go into the gas (see fig [2]). The 
velocity of a particle with equivalent velocity e z located at position z, z_(e z ) < z < 
z+(e z ), is given by 



v z (z,e z ) = sgn(e z )\l e\ - —W{z), for z G [z-(e z ), z + (e z )}, (7) 
V m 

and 

v z (z-(e z ),e z ) = v z (z-(-e z ), -e z ) =0. (8) 
Moreover, for trapped molecules we have also 

v z (z + (e z ), e z ) = v z (z + (-e z ), -e z ) = 0. (9) 

Let us define 

a z ( Zj e z ) = (el - ^(z))" 1/2 for \e z \ > y/2W(z)/m, 
so that cr z (z, e z ) v z (z,e z ) = sgn(e z ), and 

r z {e z ) = / a z (z,e z )dz= / (e 2 z - -W{z))- l l 2 dz, 

Jz-(e z ) Jz-{e z ) m 

l{e z ) = \e z \ r z (e z ). 

As in ([TTj). r z (e 2 ) can be interpreted as the time for a trapped molecule to cross the 
surface layer and l(e z ) is a length. Moreover, for every fixed z E]0,L] the application 

v z — > e z is a one-to-one application from [0, +oo[ onto [\J^W(z), +oo[ and from (JJJ) 
we have 

dv z = \e z \a z (z,e z )de z for v z > 0, e z > 0. 
On the same way, for every fixed z G]0, L\ the application v z — > e z is a one-to-one 
application from ] — oo, 0] onto ] — oo, — \J -^W(z)] and 

dv z = \e z \a z (z,e z )de z for v z < 0, e z < 0. 

Thus if we denote £ z (z) = {e z , \e z \ > (2W(z)/mY^ 2 }, then for a given function 

i/j(z,v z ) we have 

ip(z, v z ) dv z = / ip(z,v z )dv z + if)(z,v z ) dv z , 

v z Jv z >0 Jv z <0 

+ 0O 

, ip(z,v z (z,e z )) \e z \a z (z,e z )de z 



+ / ip{z,v z {z,e z )) \e z \a z {z,e z )de z 



ip(z,v z (z,e z )) \e z \a z (z,e z )de z 



ip(z,v z (z,e z )) \e z \a z (z,e z )de z 
9 



and 



/ if)(z,v z (z,e z ))\e z \a z {z,e z ) de z dz 

J£ z (z) 
z+(e z ) 

i/;(z,v z (z, e z ))\e z \a z (z,e z ) dz de z . 

e z J z-(e z ) 

To take into account the molecule-solid interaction we split the flow in two parts, 
the surface flow (for < z < L) and the bulk flow (for L < z). Since the intermolec- 
ular collisions are neglected, the state of the molecules in the bulk flow is described 
by the distribution function f(t,x,z,v x ,v z ) which satisfies the following equation 

d t f + v x d x f + v z d z f = 0. 



2.3 Kinetic model for the molecules inside the surface layer 

We consider now the population of molecules inside the surface layer and we split it 
into two groups. The group of "trapped molecules", with | e z \< ^2W m /m and the 
group of "free molecules" with | e z |> \j2W m /m. We introduce for Z-(e z ) < z < 



${t,x,z,v x ,e z ) = f(t,x,z,v x ,v z (z, e z )), for | e z \< ■ 



m 



2W m 
m 



4> f (t,x,z,v x ,e z ) = f(t,x, z,v x ,v z (z,e z )), for | e z \> J 
and 

^e z )=^ t (e z )x t (e z ) + ^(e z )x f (e z ), 

where x\ e z) is the characteristic function of the set {e z , \e z \ < ^j2W m /m) and 
X^( e z) is the characteristic function of the set {e Z} \e z \ > ^y2W m /m}. Let us notice 
that at the boundary of the surface layer [z = L) we have 

<() f (t, x, L, v x , e z ) = f(t, x, L, v x , v z (L, e z )), for e z < 0, (10) 
f(t,x,L,v x ,v z ) = (f> f (t,x,L,v x ,e z (L,v z )), for v z > 0. (11) 

Then (f> a , a = t, f satisfy 

d t <p a + v x d x <t> a --U\x)d Vx <p a + v z (z,e z )d z <p a = -(^M-f), (12) 
m r ms 7q 

where T ms is the molecule-phonon relaxation time, M is the dimensionless equilib- 
rium distribution function 
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n[(j)\ is the expression of n[f] in term of the new variables 



and 



n [0] = / / <i>(t,z,z,v x ,e z )\e g \(Tg(z,e s ) de z dv x 

'v x J£ z (z) 



70(2)=/ / M(v x ,e z )|e 2 |a 2 (z,e 2 ) de 2 dv x = exp 



»i J£z(z) 



m \ kT 



Multiplying Qnj) by \e z \a z (z, e z ), we get 



d t {\e z \az4> a ) + v x d x (\e z \a z cf> a ) ~ -U'(x)d Vx (\e z \a z <j ) a ) + e z d z <\> a = |e 2 |a 2 / p a h [0], (13) 
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where 

Let us notice that, from (EH), we have 



4>\t, x, z±(e z ), v x , e z ) = </>*(t, x, z±(e 2 ), v x , -e z ) = f(t, x, z±(e z ),v x , 0). (14) 
We introduce the following dimensionless quantities 

7 t ~ _ 7~ms ~ X _ Z _ ^ G-x/z _ V x V z . . 

t* T^s X Z V* V V 

f f 1 * ~ °* T~T U XAT W (, R \ 



where t^ s is a reference relaxation time, t* = r^ s , v* = \J (2kT)/m a reference 
velocity, x* = v*T^ ns and z* = L are reference lengths, 0* = n*/v* 2 , where n* is a 
reference number density, U* = (m/2)v* 2 , a* = 1/v*. Moreover we introduce 



n[(j)](t,x,z) = / 4>(t,x,z,v x ,e z )\e z \a z (z,e z ) de z dv x , 

'Vx J£z(z) 



(where E z {z) = {e 2 , \e z \ > yW(z)} ) so that n[<j)] = n[(f)]/n*, and 

M = exp 70 = ^ = vrexp(-iy). 

Finally we assume that 

z* 

e = — << 1, 

x* 

is a small parameter. Inserting (11511161) in equation (ITBl for a = t, f, we get 

dtde^r) +d x d x (\e z \a z 4> a ) - ^'(x)^(|e 2 |a 2 a ) + -e z d z 4> a = \e z \aj;S\, (17) 

11 



where 

T ms ^70 I 

To derive a "surface equation" we look for a solution of ( |T7|) in the following form 

where the function <p a ' t are taken so that (see ([HI)) 

^i^,x,z ± ,v x ,e z ) = 4> t '\t,x,z±,v x ,-e z ), (19) 
(p f ' l (t,x,z_,v x ,e z ) = (f) f ' l (t,x,z_,v x ,-e z ). (20) 

Inserting (|T8|) in (|T7|) . we get 

at leading order 

~eM afi = 0, 
and thus <p a, ° does not depend on z, i.e. 

^° = ^°(t,5,^,e,). (21) 

This property together with (I19H20I) implies that 

4> a '°(i, x,v x ,e z )- a '°(t, ~x, v x ,-e z ) = 0, 

i.e. implies that a '° is an even function of e z . 

at order + 1 

d~ t {\~e z \~a z ~r fi ) + v x d x (\e z \a z 4> a > ) - ^U'(x)d iix (\e z \a z 4> a ' ) + e z d~ z ^ 1 = \e z \aJ; h [<P% 

and since a, ° (and also {e^a*^' ) is an even function in e z we can take the even 
part and we obtain 

d~ t {\e z \a z fr> Q ) +v x d x (\e z \a z 4> a '°) - ^U'(x)d, x (\e z \a z ^ ) 

+ ^[d z r'\e z )-d z ~^\-~e z )} = \e z \aj; h [n 

Then we integrate with respect to z and denoting I = l(e z ) = \e z \f z (e z ) where 
f z (e z ) = J? + a z (z, e z ) dz, we get, 

<9 f (/~a)0 a '°) +v x d x (l(e z )4> a >°) - \u\x)d,Me z )4> afi ) 

12 



where 

e*[fi>](t,x,e t ) = i- [ + (n[^](i,x,z)a t {z,e t )/%(z)) dz 

T z Jz_ 

But we remark that 

*d s (F'\z,e t )-fr\z,-e z ))dz = ^([fr\z + ,e g )-fr\z + ,-e z )]- 

[0°' 1 (5_,e a )-^ 1 (i_,-e 2 )]) ) (22) 

Let us first consider the trapped molecules (i.e. a = t). From (JTHJ) we can conclude 
that this term vanishes. Thus, we obtain the following equation (in dimensionless 
form) for the trapped molecules 

We introduce the dimensionless ID distribution function g which is a density number 
of molecules per sc-unit 

g l (t,x,v x ,e z ) = l(e z )4> tfi (t,x,v x ,e z ), 

g°(t, x, v x , e z ) = /~(e z )0*'°(t, x, v x , e z )^(e z ) + l(e z )<f) f > (t, x, v x , e z )x f (e z ), 

where X*(e z ) is the characteristic function of the set {e z , \e z \ < yW^} and x^(e 2 ) 
is the characteristic function of the set {e z , \e z \ > W m }. Then the above equation 
writes 

dfg t + v x d s g t -~U'(x)d*J t = Q l ph [g\ 

where 

e[g°](t,x,e z ) = _ (n[g°/l](t,x,z)a z (z,e z )/%(z)) dz. 

Tz(e z ) J z _ 

Let us now consider the "free molecules". We first remark that (from (JTSJ) and 

dm)) 

l(e z )4> ffi = l(e z )4> f (l,e z ) + 0(e). 

Moreover the flux term is given by (122j) . with a = f and z + = 1 in dimensionless 
variables. As previously the term in S_ vanishes because of (120]) but the term in 
z + = 1 does not vanish and represents the flux of molecules between the surface 
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layer and the bulk flow (see figure [2]). 
Then 0221) for free molecules writes 



|[«p(l,g z )-<p(l,-e z )] 

which can be equivalently written (since <j>f'°(l,e z ) — 0^'°(1, — e z )=0) 

||[^' (l,e,) + e ^ 1 (l J g,) - (4> f >°(l,-e z ) +s^ 1 (l,-e z ))}. 

But from ( lT8i) . we have + ecj)^ 1 = <j>f + 0(£ 2 ), so that the flux term writes 

^f(l 7 e z )-ft(l,-e z )} + 0(e), 

and since 0-''° = $ + O(^) we have 

d- t {l{e z )ft{l,e z )) + v x d x (l(e z )^(l,e z )) - ^U'(x)d,Me z )^(l,e z )) 

+|[0/(l,g z )-0/(l,-e,)] = 



^ (0* + - ^(1, e.)) + 0(e) 



Using (TTUT) and denoting (with some abuse of notation), 

4> f (t,x,v x ,e z ) = <f> f {t,x, l,v x ,e z ), 
the equation of the free molecules finally writes 

di(l(e z )4>f) + v x d x (l(e z )$ f ) - l -tJ\3:)d,Me z )¥) 

r V(\z 



l(e z 



Introducing as previously 



+^ f (\e z \)-f(-m,e z )\)\ 



g f (t,x,v x ,e z ) = l(e z )4> f (t,x,v x ,e z ), 
g(t,x,v x ,e z ) = g\ix,v x ,e z )x\e z ) + g f (ix,v Xl e z )x f (e z ) 



this equation writes 

dfg f + v x d x g f - \u'{x)d,ji + _L[£/(|g,|) - l{e z )f{-\v z {\, ~e z )\)] = Q f ph [g] + 0(e), 
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where Q f ph [g] = (Q[g]l(e z )M - g f J. Let us remark that 0[g°] = Q[g] + O(e), so 
that the equation for trapped molecules writes 

dfg t + v x d x ~g t -^U\x)d i J t = Ql h [g] + 0{e), 

Denoting f s (e z ) = /(— \v z (l, e z )\) and neglecting 0(e) terms in the above equations 
for free and trapped molecules, we finally obtain the following system of coupled 
equations (in dimensionless form) 

dfg t + v x dfg t - l -U'{x)d,J t = QlM 

dfg f + v x d & ~g f -\u'{x)d,J f + ^^d^D -'7*1 = Q'nffl- 

Now we come back to dimensional quantities, noticing that r* = z*/v* = et^ and 
we define 

n* z* 

g(t,x,v x ,e z ) = — T cfty,x,v x , e z ), 
v 

g f (t,x,v x ,e z ) = — T g f {t,x,v X) e z ) ) 
v 

g(t,x,v x ,e z ) = g t (t,x,v x ,e z )x t (e z ) + g f (t,x,v x ,e z )x f (e z ). 
We obtain finally the system governing the flow of molecules 

Proposition 1. Under hypothesis ((TjEMP; in the limit of a small surface layer 
(e = 4r — > 0), the flow of molecules can be described by the following multi-phase 
kinetic model 

d t f + v x d x f + v z d z f = 0, (23) 

d t g t + v x d x g t --U'(x)d Vm g t = Q\ h , (24) 
m F 

WW * J U s -9 f (\e z \) _ nf 

m lr z 

l(e z )f(t,x,L,v x ,v z )\ Vz>0 = g f (t,x,v x ,e z (L,v z )), (26) 



9 tg ' + v Ag '-^, 9 '- 1 ' = Q' rh , (25) 



where 



Q a ph = — ( 



T, 



ins 



e[g](t,x,e z ) = — — / (n[g/l](t,x,z)a z (z,e z )/"/o(z))dz, 

r -Mz) J z ^(e z ) 

f s (t,x,v x ,e z ) = f(t,x,L,v x ,-\v z (L,e z )\), 
l(e z ) = \e z \ T z (e z ). 
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Remark 2. 



1. This model can be interpreted as a multiphase model. The first phase, con- 
stituting of the bulk flow of molecules outside the range of the surface forces, 
is described by a usual kinetic equation The molecules within the range 
of the surface forces are considered as a separate phase described by a two 
energy-group kinetic model, with the low energy "trapped molecules" (2~4\ ) and 
the high energy "free molecules" < f^3|) . The two groups are coupled by the col- 
lision terms. The two phases are coupled by the relation j[2h)) and the last 
term of the left-hand-side of equation If2h*\) . Let us notice that equations (Efy 
and W3\) do not give a precise description of the flow of the surface molecules 
with respect to the distance z to the surface, but give a relevant information of 
the flow parallel to the surface, and therefore will be useful for evaluating the 
transport flux in this direction. 

2. The distribution function f(t,x,z,v x ,e z ) describes the number density of gas 
molecules with velocity (v x ,v z ) in a unit (x, z)-area, but the distribution func- 
tions g l (t,x,v x ,e z ) and g*(t, x, v x , e z ) describe a number density of surface 
molecules with velocity (v x , e z ) per unit x-length. As a consequence n[0] (t, x, z 
is the number density of surface molecules per unit x-length. Moreover we can 
easily verify that 




Q P h(v x ,e z ) dv x de z = 0, 



which ensures the local conservation of mass for surface molecules. 

3. In this multiphase model, the "surface" of the solid part is identified to the 
interface z = L between the surface layer and the bulk flow. Equation j{23}) 
describes the gas flow and equations (Effi2l)\) give a simplified description of 
what can be interpreted as a motion of molecules on the "surface". The condi- 
tion If2h]) can be interpreted as a non local boundary condition for the bulk flow 
giving a description of the interaction between the gas molecules and the wall 
which is more detailed than usual local boundary conditions. More precisely, 
some molecules go from the bulk flow into the surface layer, are sent back 
by the repulsive interaction potential and immediately escape from the surface 
layer into the bulk flow, giving a specular reflexion. On the other hand, some 
molecules can go from the bulk flow into the surface layer, can have a collision 
with phonons and lose enough transverse energy so that they are trapped inside 
the boundary layer (where they are transported according to §{2J$) and, after 
some time, can gain in a collision with phonons enough energy to escape from 
the surface layer. Since r ms is small compared with the characteristic time of 
evolution of the bulk flow, this complicated interaction with the solid surface 
described by the multiphase model can be interpreted in a first approximation 
as a reflexion of the molecule by the boundary. 



16 



3 lD-mesoscopic kinetic equation for surface molecules 



3.1 Introduction 

In section 2 we derived a multiphase model describing the coupling of the bulk 
flow with the motion of molecules on the surface. This motion is given by a set 
of two coupled kinetic equations on variable x, including a Vlasov term due to the 
interaction potential parallel to the surface. In many applications this potential field 
is periodic with a small period 25 that could be much smaller than the characteristic 
distance x*, but nevertheless, for consistency with the asymptotic analysis leading 
to theorem 1, we assume that 

L « 5 « x*. (27) 

Some of the surface molecules have a small tangential energy so that they are bound 
in a potential well where they are oscillating very rapidly. At the scale x* they 
can be considered as motionless. On the contrary surface molecules with a high 
tangential energy are unbound and can move across the potential wells, but their 
velocity strongly oscillate. In so far as we are interested in the mass flow along the 
surface over a domain much longer than 5 we do not need to get information on 
the distribution function at the "microscopic" scale of a potential unit cell. Thus 
it is useful to derive a kinetic model at a "mesoscopic" scale larger than 5 but 
comparable to l mp the mean free path between two "collisions" with phonons. At 
this scale a kinetic model describing the motion of the unbound molecules by their 
average velocity over the potential cells is relevant. In this section we derive such 
a model in a simple configuration where we assume that the free molecules in the 
surface layer can be neglected so that all surface molecules are trapped. This is 
justified when W m » kT, which is true for some practical situations, and only 
very few molecules can escape from the interaction range of the normal interaction 
potential W so that + gf] « nfg*]. In such situations it is reasonable to assume 
that 

W m = +oo, (28) 

so that when a molecule enters the surface layer, it cannot escape. Moreover we 
neglect the flux of incoming molecules. This model is derived from the kinetic 
equation (j2l|) by an asymptotic analysis when 5/x* — > 0. For this asymptotic 
analysis we consider the hypothesis (jH|5]|6]), as above. Moreover we assume that the 
molecule-phonon relaxation time is constant i.e. 

T ms = const (29) 

and that U (the tangential part of the potential) is periodic with period 25. More 
precisely 

U(x) = U(x/5), (30) 
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where U(y) is a periodic potential with period 2 defined on [—1, +1] and such that 
< U(y) < U m and U(±l) = U m . 

Before to derive such a model we need to introduce some notations related to 
the motion of a molecule in a periodic potential field. We denote 



for v x ^ 0, e x = e x (x,v x ) = sgn(v x )\lv% + —U{x). 

For any fixed x the application v x — > e x is a one to one application from [0,+oo[ 
onto [^/2U(x)/m, +oo[, and also v x — > e x is a one to one application from ] — oo, 0] 
onto ] — oo, — ^2U(x)/m] so that we have 



for |e. 



> y/2U(x)/m, v x (x,e x ) = sgn(e x )\ e 2 x - —U(x). 

V m 



The jacobian of the application e x — > v x (for e x > for instance) is given by 

for v x > 0, dv x = \e x \a x (x, e x )de x , 

where 

2 

<7x(x,e x ) = (el U(x))~ 1/2 , for \e x \ > ^2U(x)/m. 

and in the same way 

for v x < 0, dv x = \e x \a x (x, e x )de x . 

Then denoting 

£ x (x) = {e X} \e x \ > (2f/(x)/m) 1 / 2 }, 
we have for every integrable function ip{ v x), 



ip(x,v x )dv x = / i/j(x,v x (x,e x ))\e x \(T x (x,e x ) de x . 
Je x {x) 

The trajectories of the molecules in the (x,v x ) plane are the level curves of the 
the total energy E(x, v x ) = me x /2 (see fig[3]). If the energy of a molecule is less than 
U m , then its trajectory is a closed curve (bound molecule), the molecule is trapped 
in a potential well. On the other hand, if its energy is larger than U m its trajectory 
is an open curve (unbound molecule) and those unbound molecules generate a flow 
in the x-direction. 

Since the tangential potential U is periodic with a small period 25, the velocity 



v x (x,e x ) = sgn(e x )J e x — -^U(x) is rapidly oscillating. If we look at the surface 
molecules on a space scale x* » 5, the average velocity of the bound molecules is 
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Figure 3: trajectories in the phase space (x,v x ) for U(y) = ky 2 /2. 



equal to zero. On the other hand the average velocity of an unbound molecule can 
be obtained as follows. We introduce —1 < y~(e x ) < y+(e x ) < +1, defined by 

U(y±(e x )) = -el for e\ < -U m , (31) 
2 m 

y±(e x ) = ±1, iove 2 x >-U m , (32) 
and we define for y E [—1, +1], 



v t{y^x) = sgn(e x )\le 2 x - —U(y), if e 2 x > ' 



m m ' 



/ 2 - 2U 

v t(y^x) = sgn(e x )\ e 2 x U(y), if e 2 x < — - and y e \y-(e x ), y+(e x )], 

V m m 

2U 

v f{y^x) = 0, if e 2 x < — - and y [y_(e x ),y + (e x )}, 

m 



and 



p ^ = sgn(e x ) 
vt (y, e x 



so that a*(y,e x ) = +oo, for y <£ [y_(e x ),y + (e x )]. Finally, we introduce a x (e x ) 
defined by 

1 ry+(ex) 1 ry+(e x ) r> 

v*(e x ) = - **(y,e x )dy = - (e 2 x - -U{y)y l ' 2 dy. 

1 Jy-M 1 Jy-(e x ) m 

Let us consider now an unbound molecule moving in the one dimensional periodic 
field U(x) = U(j) with an "equivalent velocity" e x . The position x of the molecule 
is a monotonic function of t, so that we can consider also t = t( monotonic 
function of x. Since 

dx = (e 2 x - lft£))Wdt, 
m 5 
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we have 

dt = {e 2 x - -U(%))-^dx, 

and the "time of flight", necessary for the considered molecule to cross a potential 
well (i.e. for going from x = (i — 1)5 to x = (i + 1)5) is 



Tfi(e x ) = 




where the mean velocity w x (e x ) of the molecule is given by 

. def 2 sgn(e x ) 5 

w x (e x ) = — — — . 

T f i{e x ) 

Let us notice that bound molecules, with an "equivalent velocity" e x with me 2 x < 
2U m , are trapped in a potential well, between (i + y_(e x )5) and (i + y + (e x )5). Thus 
for such a molecule, the time necessary to cross a potential well can be considered 
as infinite and by consequence the average velocity on a large scale x* » 5 can be 
considered as null, which is consistent with the above definitions of &f{y, e x ) and of 
w x . 



3.2 Homogenization of the surface kinetic model 

To describe the flow in the x-direction induced by unbound molecules we start from 
the kinetic equation (124)) for trapped surface molecules, in which we express the 
distribution functions in term of the variable e x , e z , i.e. we introduce h f given by 

h*(t, x, e x , e z ) = gf(t, x, v x (x, e x ), e z ). 

Moreover we denote x, e x , e z ) = /i*(£, x, e x , e z )x c ( e x), c — b,u, where Xb{e-x) 
is the characteristic function of the set {e x , \e x \ < ^2U m /m} and Xu(e x ) is the 
characteristic function of the set {e XJ \e x \ > yj2U m /m}. From the definition of 
y±(e x ) ( 13TH32|) . it comes for bound molecules 

hl(t,Sy ± (e x ),e x ,e z ) = h\{t, 6y±(-e x ), -e x , e z ) = g\t, 5y±(e x ), 0, e z ). (33) 

Let us recall that the evolution of the surface molecules is described by equations 
( l24l) which can be written 

d t g t + v x d x g t --U\x)d v J = Qt h , 
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From the definition of h l , we get d t h l = d t g t ', and d x Y^ = d x g l + d v g t d x v x , or 

mv x (x,e x ) 



d x h l = d x g l u ^ s d^q 1 , so that 



v^x.e^dxh 1 = v x (x,e x )d x g t - —U'(x)d Vx g t , 

and finally the distribution function h l {t, x, e x , e z ) satisfies the following kinetic equa- 
tion 

dth' + v^x.e^h 1 = — (9[tf]Z(e z )M - h l ) , (34) 

Tins 



where 



Q[h] = / (n[h/l}a z (z,e z )/'ji(x,z))dz, 

Tz{e Z ) J z _(e z ) 



n 



/ / v(t,x,e x ,e z )\e x \a x (x, e x )\e z \a z (z, e z ) de x de z 

J £ z (z) J 8 x (x) 



,{z) J £ x (x) 
-m( e 2+ e 2)/(2fcT) 



M(e x ,e z ) = e 

, , 2nkT ( U(x) + W(z) 
7i(x,z) = exp 



m \ kT 

To derive a kinetic model at a "mesoscopic scale" (i.e. on a characteristic length 
x* » 6), we introduce the following dimensionless quantities 

~ X ~ I _ V x ^ G-x/z ~ t _, T ms „ 77- /„ p.\ 

37 — . I — . U x — . &x/z — i ^ — 5 ^~ — ; ^ — ; V J 

x* z* v* v* t* t* n* 

and 



where x* is a reference length, t>* = ^2kT/m, t* = x*/v* is a reference time, 
U* = (m/2)v* 2 , n* is a reference number density, h* = (n*z*)/v* 2 , a* = 1/V, so 
that 

v x (x,e x ) = sgn{e x )\Je 2 x -U(xx*/6). 
Moreover we assume that the reference length x* » 5, so that the ratio 

S 

e = — << 1 

x* 

is a small parameter. We now derive an asymptotic model in the limit e — > 0, and 
we describe the asymptotic expansion for the trapped molecules. 



Inserting (I35H36I) in the kinetic equation (134|) we obtain 

<%(/?) + sgn{e x )\le 2 x -U{x/e)d x {h t ) = - (e[#]ZM - h* 
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and we look for a solution in the following form 

h^t, x, -, e x , e z ) = h tfi (i, x, -,e x , e z ) + eh 1 ' 1 ^, x, -, e x , e z ) + ... (37) 

where h t,l (t, x,y,e x ,e z ), i = 0,1,... are periodic functions in y with period 2 and 
where, for bound molecules (to be consistent with f[3"3"j) and the definition of y-(e x ) 
and y+(e x )) 

hb\t, x, y±(e x ), e x , e z ) = h\\t, x, y±(-e x ), -e x , e z ). (38) 
Then h l satisfies 



dffi) + sgn{e x )y]<? x -U{x/e) [d x + U v \(h*) = \ (e[h*\lM - ti) , (39) 

where 

e[#] = \- [ +{ ' {n\ti/l\ a z (z,e z )/-yi (x/e,z))dz, 

T z Jz-{e z ) 

nW\ = / ^\t,x,y,e x ,e z )\e x \af(y,e x )\e z \a z (z,e z ) de x de z . 

J£ z (z) J£#(y) 



£ z (z) = {e z ,\e z \>JW(z)} 



£f(y) = {ex,\e x \ > \/U{y)} 
7i (y, z) = 7rexp \-U{y) - W(z)j . 

Moreover since 



v x {x,e x ) = sgn{e x )\Je 2 x - U(x/e) = v*{^,e x ), (40) 

inserting (13 7} |4"01) in f)39p and balancing order by order in e we get at the principal 
order, 

Vfdyh^ = 0, 

and thus for molecules with vf ^ we conclude that 

/i*'°(t, x, y, e x , e z ) = h tfi (t, x, e x , e z ). (41) 

Then from this property and from (!38|) we deduce that h% (t, x, e x , e z ) is an even 
function of e x . 

At the next order we get 

dfh 1 ' + vfd^ + vfdyti* 1 = i (&[h°}lM - . (42) 
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We consider now separately bound and unbound molecules. 

Let us first consider unbound molecules. Taking into account ()4ip and multiply- 
ing (J4"2j) by e x ) _1 , we get, for unbound molecules 

— dfh^ + ajtf + d y h^ = 1 (e[^°]/M - I'A . 

v£{x/e,e x ) v£(x/e,e x ) r \ ' 

Then averaging with respect to the fast variable over one period and taking into 
account that h^ 1 is periodic in y, we get 

1 ^.r, 



, dirt + dirt = — — -= (eih^m - h$°) , 

w x {e x ) w x (e x ) t V 



where 



w x {e x ) = sgn(e x )(^j v*(y,e x ) 1 dy^j 



1 



2a x (e x ) 



+i 



= / ^(y,e x )e[h^](t,x,y,e z ) dy. 



i 



Finally, multiplying by w x (e x ) we obtain 

afif + «J.(^)flS6^° = i (pp°]iM - ~rt°) , (43) 

Let us consider now bound molecules. Since h^°(t,x,e x ,e z ) is an even function 
of e x it is equal to its even part and thus the left-hand-side of (142]) writes also for 
bound molecules 

+ \(vt(e x ) + vf{-e x ))djif + \ (v*(e x )d y h^(e x ) + ^(-e^dyh^-e^) . 

But vf(y,e x ) is an odd function of e x , so that (vf(e x ) + vf(—e x ))d x h^° cancels. 
Then multiplying by af(x/e, e x ) and averaging in the fast variable over one period 
(between y-{e x ) and y + (e x )), we get finally for the bound molecules 

dthi° = \ (m tfi }lM - hi ) , (44) 



where 



— ~ 1 rv+i^x) 

0[^'°] = —7-, T^rr^T / affaejeih^foxMe,) dy. 



(y+(e x ) - y-(e x ))cr x (e x ) J y _(e x ) 
Rewriting equation (I43H44I) in dimension form and omitting the superscript 0, 



we get 



d t h t u (t,x,e x ,e z ) +10^)8^1^, x,e x ,e z ) = — (G[h*\lM - ft) , (45) 

d t hi(t,x,e x ,e z ) = — (Q[h*\lM - h$ . (46) 
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Proposition 2. Under the assumptions [$M M2l§28\\29fi3(%) . the solutions of ffift) 
formally converge as e = 5 / x* — >■ to a distibution function ^(t, x, e x , e z ), satisfying 
the following "mesoscopic" kinetic model 



dtti + w x (e x )d x h t = J- (e[h*\l(e s )M - h*) 



where 



w x {e x ) = Xu{e x ) sgn(e x ) - / (e 2 x - 2U(y)/m) 1/2 dy 



2 



i 



&[tt\{e x ,e z ) = - — j-^r / a*{y,e x )Q[h t ](y,e z ) dy, 

(y+(e x ) -y-{e x ))a x {e x ) J y _ (ex) 

@[^](y,e 2 ) = —— / (n[h/l](y,z)a z (z,e z )/^i(y,z))dz, 

, , 2nkT ( U(y)+W(z) 
1i{y,z) = exp 



m \ kT 



n[g](y,z) = / g(e x , e z )\e x \a*(y, e x )\e z \a z (z, e z ) de x de z 

JEJz) Jet(y) 



>S z (z) J£?(y) 

*x{ex) = o/ \\ e l — u (y)) d y> 



2 Jy-(e x ) V V m 
S*(y) = {e x , e 2 x > 2U(y)/m}. 

Remark 3. 1. The distribution functions h\, for c = b,u describe the number 
density of gas molecules with velocity (e x , e z ) per x-unit, obtained by averaging 
the distribution functions g t over the periods of the tangential potential U. 
Then 

ry-i e =o) r-L p p 

N(t,x)= / / n[h/l}(t,x,y,z) dz dy = 2 / h\e x \a x (e x )de x de z 



Jy-{e x ) JO Je z J e x 

is the number density of molecules per unit-length, averaged over the periods 
of the tangential potential U. 

2. We can check that J e J e (Q[h]lM — h)\e x \a x (e x )de x de z = 0, which ensures 
the local conservation of mass. 

3. In this model the unbound surface molecules move in the x-direction at an 
average velocity depending on their total tangential energy, or equivalently, on 
e x . On the other hand the bound surface molecules are trapped in a well of the 
tangential potential U. Their average velocity is null and they relax inside a 
potential well toward the equilibrium (see equation Let us notice that, 
contrary to ([TBj). we do not need to assume that the bound molecules are at 
equilibrium. 
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4- We have assumed above that L « 5 « x* (21). It could be interesting to 
drop this hypothesis and to consider the case where L and 5 are comparable. 
But, in such a case, we cannot uncouple the asymptotic analysis of propositions 
1 and 2 and we have to treat both of them together. 



4 Diffusion models for surface molecules 

In section 2 we derived a multiphase model describing the coupling of the bulk flow 
with the motion of molecules on the surface. This motion is given by a set of two 
coupled kinetic equations on variable x, including a Vlasov term due to the inter- 
action potential parallel to the surface. In this section we will consider integration 
time much greater than the molecule-substrate relaxation time and thus we derive 
diffusion models. We consider the case of a smooth interaction potential parallel to 
the surface. We proceed in several steps and we begin with simplifying assumptions. 

In a first subsection we consider the same simple configuration as in section 3 
where we assumed that the free molecules in the surface layer can be neglected so 
that all surface molecules are trapped. As before it is reasonable to assume that 

W m = +oc, (47) 

so that when a molecule enters the surface layer, it cannot escape and moreover 
we neglect the flux of incoming molecules. Thus there are no free molecules inside 
the surface layer and it is possible to describe the trapped surface molecules by 
the following closed model describing the evolution of the distribution function g = 

g(?x, e 2 ) 

d t g + v x d x g - ^d Vx g = — (Q[g]lM - g) , (48) 

TTl Tins 



where 



and we denote 



M{v x ,e z ) = ex P (-m(v 2 x + e 2 z )/2kT) 



mvl \ , , , s / met 



M x (v x ) = exp -— £ , M z (e z ) = exp 



2kT 7 K ' L \ 2kT 



2knT . . f %r . ,. , , , 2knT ( W(z) 
lx = \ , lz{z) = / M z (e z )\e z \a z {z,e z )de z = \ exp 



m J Sz{z) V m \ kT 

The diffusion limit is derived, first in the isothermal case and then extended to the 
non-isothermal case. 

Afterwards, in the second subsection we drop the assumption W m = +oo and we 
extend the analysis of the diffusion limit to the configuration of a flow in a narrow 
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channel where the free molecules are not neglected but where there is no bulk flow. 



Before deriving the diffusion models we need a technical lemma related to the 
collision operator Q[g]lM — g. Let us denote in the following, 

((•)) = / dv x de z . 

Je z Jv x 

Then we have 

Lemma 1. Let us consider the following equation for g: 

e[g](x,e z )l(e z )M(v x ,e z ) - g(x,v x ,e z ) + ip(x,v x )l(e z )M z (e z ) = 0, (49) 
where ip is a given function. 

A necessary and sufficient solvability condition for this equation is 

ij)(x, v x ) dv x = 0. (50) 

Moreover if /i50\) is satisfied every solution of fijty ) writes 

g = /3(x)l(e z )M(v x , e z ) + ip(x, v x )l(e z )M z (e z ), 
where (3 does not depend on v x and e z 

Proof, (i) As we have noticed in remarks 1, ((Q[g]lM — g)) = which gives the 
solvabiliy condition fl50|) . 

(ii) It can be easily checked that if /3 = /3(x) does not depend on v x and e z , then 
Q[f3lM] = (3. Thus g(x,v x ,e z ) = /3(x)l(e z )M(v x ,e z ) are solutions of equation ( 149]) 
with if; — 0. We can prove that they are the only solutions. Indeed if we look for 
a solution in the form j3(x,v x ,e z )l(e z )M(v x ,e z ), then /3(x,v x ,e z ) is solution of an 
integral equation 

/3(x,v x ,e z )= / k(v Xj e Zj v' Xl e' z ) /3(x,v' xj e' z ) dv' x de' zJ (51) 

Je' z Jv' x 

where the kernel k(v x ,e z ,v' x ,e' z ) = k(e z ,v' x ,e' z ) is defined by 

, ,s r +(e *' ei) <T z (z,e z )\e' z \a z (z,e')MJv' x )M z {e') , 
k(e z ,v' x ,e' z ) = / V n ' \\ z) *\ J V zl dz, 

Jz-{e z ,e' z ) lx lz\z) r z {e z ) 

where 

z_(e z ,e' z ) = m&x(z_(e z ),z_(e' z )), z + (e z ,e' z ) = mm(z + (e z ), z+(e' x )), 
and this kernel k satisfies 

V e z , / / k(e z ,v' x ,e' z ) dv' x de' z = 1. 
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Moreover, since the kernel k is independent of v x , so is /3 and the integral equation 
(I5T1) writes 

(3(x,e z )= [ k(e z ,e' z ) 0{x,e' z ) de' z , (52) 



where the kernel k(e z ,e' z ) = J v , k(e z ,v' x ,e' z )dv' x satisfies 

V e z , k(e z ,.) e L\, z and / k(e z ,e' z ) de' z = 1. 

Je' z 

It is easy to check that functions = 0(x) (constant with respect to e z ) are solutions 
of this equation. Moreover we prove now that they are the only solutions in L^. First 
we prove by contradiction that the only solutions of a(e z ) = j e , k(e z ,e' z )a(e z )de' z 
are constant functions. If we denote m = ess sup e /_a(e^.), since a is not a constant 
function of e' z 

3r] > 0, a(e' z ) < m — r), a.e. in e' z on a set JF X)TI with positive measure in e' z , 
so that 

a ( e z) = / k(e z ,e' z )a(e' z )de' z , 



< m-rj J k(e z ,e z )de z <m-r]k , 



where 

This last inequality is in contradiction with the definition of m so that we conclude 
that necessarily a does not depends on e z . We easily deduce of this property that 
the only solutions of (lo"2"|) are functions 0(x, .) which, for almost x, are constant 
functions in e z . Indeed, for £ given in C^, let us denote a^(e 2 ) = f x £(x)0(x,e z )dx. 
Then is solution of the integral equation a^(e z ) = f , k(e z , e' z )a^(e' z )de' z , so that 
is a constant. But £ — >■ is a linear form that can be written J ((x)£(x)dx, 
so that almost everywhere in x we have 0(x, e z ) = ((x). Thus every solution of the 
homogeneous equation fT^9|) writes 0(x)l(e z )M(v x ,e z ). 

(iii) Let us now consider the nonhomogeneous equation. First, we remark that 
floUl) implies that 



^ ^(x,^)^^ ^ |e z |cr z (2;, e z )M z (e z )de z j = 0, 



so that ^>£M is a particular solution of (149]) . and then, from (ii), every solution writes 

g(x, v x , e z ) = 0(x)l(e z )M(v x , e z ) + ip(x, ^(e^M^e*). 

□ 
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4.1 Diffusion model for the trapped surface-molecules 

For the sake of simplicity, we assume in a first step that the temperature T is a fixed 
given constant. 

4.1.1 The isothermal case 

Under the above assumptions the kinetic model obtained in the previous sections 
writes 

d t g + v x 8 x g - ^-d Vx g = — (Q[g]lM - g) . (53) 

To obtain a dimensionless form of this equation we introduce the following reference 
quantities, z* = L is a reference length in the z-direction, v* = ^2kT/m is a 
reference velocity, t* d is a reference diffusion time, t^ s is a reference molecule-phonon 
relaxation time, t* = z*/v* is a reference time for crossing the surface layer, n* 
is a reference number density and g* = (n*z*)/v* 2 , U* = mv* 2 and M(v x ,e z ) = 
exp (— (v 2 + e 2 )). Moreover we introduce another reference time t* such that 

< « r* ms « t* « t% 

and the small parameter e, and we set 

* ±* 

^ ms _c_ 

t* t* ' 

Finally we introduce reference length in the x-direction defined by 

x** = v*t*. 
We rescale the problem according to 

~ % ~ Z ~ ^ _ ^~rns ~ Vx ~ _ &x/z ftz.A\ 

OC — . <2> — . t — . Tj — . Tyxis — 5 ^ x — 5 ^x I z — 5 I " / 

a;** t^ s v* v* 

and 

- 9v* 2 (J= U_ q = ^1 m 
n z U n 

With the previous notations and the above assumption, the dimensionless form of 
equation fl53l writes 

edfg + v x d x ~g-U\x)dvJ = J- (e[g]lM - g) , (56) 

We consider the asymptotic analysis of this problem when e tends to 0. We look for 

~g = ~g {0) + e~g^+e 2 gV + (57) 
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Inserting ( 15 7p in ( 15 6 p and balancing order by order in e we get 
At leading order 

g(o) = e[^(°)]ZM. (58) 
From the dimensionless form of Lemma 1 the solutions of (l58~l) write 

g {0) = a(t, x)lM = N0(t ' x) lM, (59) 

7 

where N° (t, x) is defined by 

N°(t,x) = a(t,x) 7, 
and where the constant 7 (depending only on W and T) is defined by 

7=/ / M(v x ,e z )l(e z ) dv x de z . 

J e z Jv x 

Then we get 
at order +1 

v~ x d x ~g^-U\x)dv x ~g^ = J_ (e\gW]lM - fif«) . 

We notice that the left-hand-side of this relation writes 

v~ x d x N^4^- + 2U'(x)v x NM x (v x ) ) l(e z )M z (e z ), 

so that we can use the dimensionless form of Lemma 1 and thus solutions of this 
relation write 



= a 1 ft, x)lM - ^ f d x N°(t, x)) lv x M + I^if/'(x)(^iV (t, x))lM. (60) 
7 V / 7 

g w = a\t, x)lM - ^ (d x N°(t } x)) (v* - U'{x))lM. (61) 



at order 



dfg^+ v x d x ~g^-U'(x)dv x ~gV = J_ (e[S< 2 >]ZM - g®) . 

As above and taking into account (|6"TT) . the left hand side of this relation can be 
written in the form ^(x,v x )l(e z )M z (e z ), and thus existence of a solution in g^> is 
ensured under the solvability condition 

dfg m + v~ x d x g^ - U'(x)dv x g^ de z dv x = 0. (62) 
29 



Inserting ()59H6ip into (I6"2j) and noting that J\. f- U'(x)dv x g^de z dv x = 0, we get 

d;N° - D^dj 2 N° - f ms d x (u'(x)N ) = 0, (63) 
where the diffusion coefficient Dq is given by 



D (n) _ Tms_ / / ^2 M(v x ,e z )l(e z ) dv x de z . 

7 Je z 

Then we come back to dimension quantities and we denote N(t,x) = n*z*N°(t,x). 
With this definition g° = (Nfi)lM, where 7 = ((M)), and 



((.}) = / -dv x de 

J e z J e x 

so that the number density of molecules 

77 -[5 , °/^] = / / ,/ r 2 |e 2 |cr z (2:, e z ) d^rfe^ 

•/£*(*) Jv x l ( e z) 

satisfies 

y nfo/Z] cfe = / / M{v x ,e z )dv x de z = N{t,x). 



Z " U X 



Thus N(t, x) can be interpreted as the number density of molecules per unit x-length 
in the surface layer, obtained by integrating (with respect to z) over the width of 
the layer and from (I6"3"|) we get 



d t N - D^dl 2 N - r ms d x ( ) = 0, 



U'{x) 
m 



Finally we have obtained the following result 



Proposition 3. Under the hypothesis ((TjEMHlSjT^ )> ^ e solutions of Hlfy formally 
converge as e(= t*/t* d = r^ s /t*) — > to (N(t,x)/j) l(e z )M(v x ,e z ) where the function 
N(t, x) is solution of the following diffusion equation 



d t N - D^ ] d 2 ,N - r m .A l ^^N ) = 







U'{x) 



m 



where 

= ^((v x 2 IM)), 7 = «M», M(v x ,e z ) = eM*+4)/w. 
7 

Remark 4. 1. This diffusion equation describes the gas flow in the limit e — > 
only on a long time scale « td, or, which is in some sense equivalent, for 
an initial distribution of the form (N/^y)M(v x ,e z ). It does not give a correct 
description of the flow for a general initial condition and on a short time scale. 
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4.1.2 The non-isothermal case 



In this section we assume that the temperature is a (given) function of x, and 
we perform the same asymptotic analysis as in the previous section. We only in- 
dicate the modification induced by this assumption. We first notice that M = 
exp (—(v 2 + e 2 z )/ (2kT(x)) and that 7 = j(T(x)). Moreover the derivative of 7 with 
respect to the temperature, denoted by 7', is given by 

l'(T) = d Tl = ^ f - 2 (((vl + el)M}}. 

The asymptotic analysis follows the same steps as in the previous section. The 
first difference occurs at the identification of the expansion at order +1. Taking into 
account the variation of T with respect to x, the first term in the expansion reads 
now instead of (|61[) 

+f ms U'(x)d ij J , 

where the effect of the temperature gradient is included. Finally the asymptotic 
analysis leads to the diffusion equation 

d t N - d x (D^d x N) - d x (D<pd x T) - r ms d x (^-Nj = 0, (64) 

where 

D (n) = D^\t) = -^((v 2 JM)), 

V ' ; 7 (T(x)) Ux V 2kT(x) 2 l 2 {T{x))J " 

Of course, this approach can be easily extended to the case where the tempera- 
ture field is not given but determined by a conduction equation (for instance). 

It is sometimes preferred to write the diffusion equation (164]) with the mass flux 
defined with respect to d x p rather than with respect to d x N. Assuming that the 
pressure of the gas is given by the ideal gas law p = kNT, the diffusion equation 
( l64"l) writes 

d t N - d x (Ci p) d xP ) - d x (Cjpd x T) = 0, 

where 



kT 
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Figure 4: narrow channel with two surface layers. 



4.2 Diffusion model for trapped and free molecules in a nar- 
row channel 

In the previous subsection we derived a diffusion limit of the tangential model under 
the assumption that all molecules in the surface layer are trapped. Now we drop 
this limitation so that the flow of trapped molecules is coupled with the flow outside 
the surface layer by the free molecules. We consider the following configuration: the 
gas molecules move in a narrow channel with diameter 2L, with a lower boundary 
(denoted 1), located at z = and an upper boundary (denoted 2), located at z = 2L 
(see figure H]). In this simple configuration there is no bulk flow and the incoming 
free molecules in the lower surface layer are the free molecules going out of the 
upper surface layer and conversely. We assume moreover that the two boundaries 
are similar (same temperature, same interaction potential) so that 

T 1 = T 2 , U 1 (x) = U 2 (x) = U{x), W 1 (z) = W 2 (2L-z). (65) 

Using the same approach as in section 2, the flow of molecules in the channel 
can be described (in the isothermal case) by the following system of two coupled 
monodimensional kinetic equations 

d t 9i + v x d x g x - ^-^-d Vx gi = Q ph [gi] + * (g 2 (-\e z \) - gi(\e z \)) , (66) 
m 2r z (e z ) 

d t g 2 + v x d x g 2 - ^—^d Vx g 2 = Q ph [gx] + * (gx(\e z \) - g 2 {-\e z \)) , (67) 
m 2r z {e z ) 

where g\ is the distribution function describing the (trapped and free) molecules 
inside the lower surface layer and g 2 is the distribution function describing the 
(trapped and free) molecules inside the upper surface layer. From hypothesis (165]) 
the times t z is the same in the two equations and we have the same property for 
T ms , for the length /, for the characteristic function x , f° r the operator O and the 
distribution M. 

We denote g* — gi + g 2 , the distribution function of molecules in the channel. 
Then if we add equation fl66|) and equation (|6"7]) . the sum of the coupling terms 
vanishes since the free molecules outgoing from the lower surface layer are the free 
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molecules incoming into the upper surface layer and conversely. Thus (using that 
the operator is linear in g), we obtain for the following closed kinetic equation 

U'(x) 1 
d t g* + v x d x g* d Vx g* = Q ph [g*\ = (Q[g*]lM - g*) , 

which is very similar to (153j) . but includes the free molecules (gi = x t 9i + 9% , i — 
1,2). If we rescale this equation as in (l54H55p . we get the following dimensionless 
equation for g* 

+ v x d x g* - U'(x)dv x g* = — (®[g*\lM - g, 

and the same asymptotic analysis as in section 4.1 leads to the following diffusion 
equation for iV* = J n[(gi + g 2 )/l]dz 

d t N* - D^N. - r ms d x (^-N*) = 0, (68) 

where 

Dt l) = —((v x 2 lM)), 7 = ((M», M(v x ,e z ) = e -M*+*)/™ . 
7 



Let us look now for a diffusion limit of equations ( 1661I671) . We need some more 
assumptions. We rescale these equations as in ( 1541551) with, in addition 



~ _ Tz_ ~ f ~ X f 9i 
T* (rif/v* z ) 

where r* = z*/v*. We denote 

t* _ z* 
t* x* ' 

and as in 4.1.1 we assume that 

r* tr* 

t* t* 

With these notations (1661) writes in dimensionless form 

edigi + v x dr x ~g x - U\x)d^~g x = - (e[gx]lM - g^j 

and we have a similar formula for ( |67|) . The right-hand-side of those two equations 
contains a collision term (of order i) and a coupling term which makes the two 
equations relax one towards the other. To go further we must precise the relative 
size of the coupling term with respect to the collision term which is determined by 
the ratio ^-4 = Different regimes can be encountered according to the size of 
n*j:/n* compared with the small quantities e and Eq, and we study several ones in 
the following. 
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4.2.1 Strong coupling of the two surface layers 

We assume there that 

. n* f e 
n e 

The first assumption means that r* < r^ s . The second assumption means that 
the ratio of free molecules is of order — , so that the number of free molecules can 
be (smaller than but) comparable to the total number of molecules (when e Q = e). 
This is a reasonable assumption when W m ~ kT. Under those assumptions the 

ratio ^f-| = l/£ so that the coupling term of the two equations is strong and is 
comparable to the collision term. Then the dimensionless form of the system is 
given by 

efhgi+Vxdsfc-U'ffidijj! = A[gi,h}, (69) 
£dtg2 + v x d £ g 2 -U'(x)d Bm g2 = A[<? 2 ,<7i], (70) 

where A[g a ,g b ] = -^k- (&[g a }lM - g a ^ + (g b (-\e x \) - g a (\e z \)). 

Since we have derived a diffusion limit for g* = g\ + g 2 , in order to get the diffusion 

limit of the system of equations modeling the evolution of g\ and g 2 it is sufficient 

to study the diffusion limit of the equation giving the evolution of g** = g\ — g 2 

obtained by subtracting (1701 from ( 1691) . Looking for g^ = g^ + eg\^ + inserting 

in (1701) -(1691) and balancing order by order we obtain 

At the leading order 

e[g°jM-~gl) - ^j^(gl(\e z \) + &(-|S,|)) = o. 
J I r z {e z ) 

We first remark that this relation implies that g^ is an even function of e z , so that 
it can be written 

Q[g°jM - g°J - fm : X !^ z) ~gl = 0. (71) 

To solve this equation we proceed as in the proof of Lemma 1 and we look for 
g** = /3(x,v x ,e z )lM. Inserting in ( 1711) we find for /3 the following integral equation 

1 + ^yjjy^) P&VxiZ s ) = f [ ^v XJ e z ,w XJ e z )(3(x,w x ,e z )dw x di z , (72) 



€ z J w x 



where the kernel k ( the dimensionless form of the kernel k introduced in Lemma 1) 
is nonnegative and satisfies 

k(v x , e z , w x , e z )dw x de z = 1. (73) 




Thus from (172"]) . we deduce that (3 vanishes or cannot be independent of v x and e z 
and moreover 

min/3(x,w x ,e 2 ) < ( 1 + ^^0^] (3(x,v x ,e z ) < m&x/3(x,w x ,e z ). 

w x ,e z y T z {e z ) J w x ,e z 
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Since ^1 + Tm ~ x ^S e ^ j > 1 for \e z \ > W m , this relation implies that /3 reaches its 
maximum for a value (w*, e*) such that |e*| < W" m . Then we have 

/3(x,C*,e*) = / / k{vl,el,w x ,e z )(3{x,w x ,e z )dw x de z . 

J € z J W x 

But this relation and ( 173]) imply that /3 is independent of u> x and e z , and thus, from 
the remark above, that (3 = 0. Consequently we have proved that = 0. 
Using this result we check that g\ m satisfies the same equation as so that we 
can also conclude that g\^ = and so on. Finally, since g^ = 0, we conclude that 
N\ = N 2 = N*/2 and that Ni and N 2 satisfy the same equation as iV* f )68|) . In fact 
since the coupling is strong, the distribution functions g± and g 2 relax very fast to 
the same limit N if /{2^f)lM , in a time shorter than the diffusion time. 



4.2.2 Moderate coupling of the two surface layers 

We assume now that 

n *f 
n* 

n*,t* 

so that the ratio -£-4 = 1. This assumption means that the number of free molecules 
is much smaller than the total number of molecules, which is reasonable when 
W m >> kT. Under this assumptions the coupling term is moderate (of order 1) and 
smaller the collision term (of order 1/e). 
Then the dimensionless form of the system is given by 

edigx + v x d x g! - U'^d^gx = A[g l7 g 2 ], 
ed- t g 2 + v x d x g 2 - U'{x)d^g 2 = k[g 2 , g ± ], 

where 

^\3a, h] = (®[9a]lM -~g a )+^r (g b (-\e z \) - U%\)) ■ 

We look for g { = g^ + eg] + with N[cji} = n[gf/l]dz = 0, for i = 1,2 and k > 1. 
The asymptotic analysis leads to 

At the leading order For i = 1,2, Q[gf\lM — g® = 0, which implies = ^IM. 
At order +1 

2r z {e z ) V J 

A necessary condition of solvability is that the integral of the right- hand- side with 
respect to v x and e x vanishes, which implies that 

(N 2 -Ni) f [ ^^-l(e z )Mdv x de z = 0, 
Je z Jv x 2r 2 (e z ) 
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and thus N\ = N 2 . With this condition the hypothesis of Lemma 1 is satisfied 
by the right-hand-side and taking into account that N[gl] = , we obtain g\ = 
v x dx(ji ~ U' (x)d% x gf . Finally, coming back to dimension variables we conclude that 
Ni = N 2 = N*/2 are solution of a the same diffusion equation as iV*. 



4.2.3 Weak coupling between the two surface layer 

We assume now that 

Tl* 

4 = ££ o- (74) 
n* 

This assumption means that the ratio of free molecules is very small. Under this 

n*t* 

hypothesis the ratio — f4 = £, so that the coupling of the two equations is weak 
(the coupling term is of order e while the collision term is of order 1/e) and the 
dimensionless form of the system governing the evolution of gi and g 2 is 

edigx + vAgx-U'^d^Jx = A[gi,g 2 ], (75) 

£d i g 2 + Vxdxg2-U'(x)duJ 2 = A[<7 2 ,<7i], (76) 

where 

M9a, U = -4- (e[g a )lM - ~g a ) + (ft(-|e,|) - &(|e,|)) . 
£ T ms V / 2T z (e z ) 

To study the diffusion limit of this system we integrate the two equations with 
respect to e z . It comes 

edfh + Vxdih - U'fflds.fo = — ^- I ^M x -hA+Uh- Fi) , (77) 

£ T ms \ lx ) 2 V / 



edfh 2 + v x dx~h 2 - U'(x)d,h 2 = — L_ [ ^M x -h a )+UP 1 - F 2 ) , (7* 



where hi = f iz 9ide z , N { = N[gi] = N[hi], Fi = ^f L g i {e z )de z . Unfortunately, 

this system is not closed because, in general, we cannot write Fi as a function of hi. 
Nevertheless 

hi ~ ~ ~ Je X s \e z )\e z \M z (e z )de z 

if gi = — — — —l(e z )M z (e z ), then F { = c hi, where c 



/ lM z de z J lM z de z 

(79) 

But an asymptotic analysis from the dimensionless form of ( 17511761) proves that g^ = 



(Ni/y)lM and ~g\ = {Ni/tylM - ^ [d s N? vJM - U\x)d,m, so that $ + eg} 
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satisfies the assumption of (1791) . Thus we can write (I77H78P in the following form 
edfh + v x djn - U'^dvM = -L- I ^M x - h ) + (h 2 - h) + 0(e 2 ) 

S T mx \ It 2 \ J 



£ Tms \ lx 
£7~ms \ lx 



£dfh 2 + v x dzh 2 -U'(x)dz x h 2 = ^I^M.-h^+^dfh-h^+Oie 2 ) 

er m .x \ It 2 \ J 



where j x = J- M x (y x )dv x = y/2n. The asymptotic analysis on this system follows 
the same ideas as in previous subsections and we get 
At the leading order 

A/ -0 ~ _ „ _ 

h° = -*-M x , 2 = 1,2, with jVf + N° = N°. (80) 

% 

At order +1 



N 1 ~ N° 
h\ = ^M x -f ms \d x -^v x M x -U'(x)d, x h° i 

lx \ lx 

At order +2 

dfhl + v x d x hl-U'(x)d,M = J-i^Mx-h^+Uhl-h* 



lx 

dfhi+v x dihi-u'(x)d, x h l 2 = ^-[?lM l 4|]+5(i;-ii 



With this form of the equations, we can use Lemma 1 and h\ and h 2 , can be defined 
from those relations if and only if the following solvability conditions holds 



After inserting (180118 1J) in the above relations and coming back in dimension variables 
we finally obtain the following result 



Proposition 4. Under the hypothesis (RHMMlMSSf^p, the solutions g^ i = 1,2 of 



{6B§6l\ ) formally converge as e{— t* c /t* d = T^ s /t*) — > to 

(Ni(t,x)/j)l(e z ) M(v x ,e t ),i = l,2 
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where the functions Ni(t,x),i = 1,2 are solutions of the following system of diffusion 
equations 



d t N 2 - D^d 2 x ,N 2 




c {N 2 -N x ), 



c (tfi-tfa), 



where 



o 



ma 



((v x 2 M)), 



7 



T = «M», M(^,e 2 ) 



e -m(^+^)/2fcr j c = c{Wm) 



5 ez X 1 ' (e z )\e z \M z (e z )de z 
J e J(e z )M z (e z )de z 



Remark 5. From this model, thanks to convenient rescalings, we can recover some 
of the diffusion models obtained above 

1. In the limit of a large W m the number of free molecules tends to and so does 



the coefficient c = c(W m ) = (J ez X f (e z )\e z \M z (e z )de z )/(J e J(e z )M z (e z )de z ). 



Thus we recover in this limit the result of proposition 3, where we found a 
diffusion equation for the trapped molecules assuming that the free molecules 
could be neglected. In the present configuration, we obtain, in the limit of a 
large W m two independent diffusion equations for N\ and N 2 . 

2. In the limit of a small t z , then N 2 — N\, and we recover the diffusion model 
of the moderate coupling regime. 

5 Conclusion 

We have presented the formal derivation of a hierarchy of models describing a gas 
flow in the vicinity of a wall, using scaling and systematic asymptotic analysis. Fol- 
lowing the ideas introduced in [7],[8],[1], [15] , [5] , [16] , [6] , the influence of the wall is 
taken into account through Van der Wall forces acting on the gas molecules and 
through a relaxation of the gas molecules by the substrate. 

In this paper we made some simplifying assumptions: we assumed that the 
molecules move in a 2D plane, we considered the case where the intermolecular 
collisions can be neglected and we assumed that the interaction potential has a sim- 
plified structure. 

With those assumptions we derived a multiphase model involving a classical kinetic 
equation for the bulk flow coupled with two one-dimensional kinetic equations mod- 
eling the trapped and free molecules inside the surface layer. This one-dimensional 
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kinetic model can be interpreted as a non-local boundary condition for the bulk 
flow. 

Then, assuming that the interaction potential is rapidly oscillating in the direction 
parallel to the solid surface, an averaged mesoscopic kinetic model is obtained by 
homogeneization. 

Finally, in the limit of a small relaxation time, we derived from the multiphase ki- 
netic model diffusion models for the surface molecules. In a first step we assume 
that the free molecules can be neglected and in a second step we consider a narrow 
channel constituted by two surface layers. Then we took into account the trapped 
and free molecules in the channel and we derived several diffusion models according 
to the (strong, moderate or weak) coupling of the two surface layers. 
The extension of those models for more general interaction potentials and for colli- 
sional flows will be studied in forthcoming papers. 
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